Biophysics and Physicobiology

To understand the regulation of gene expression, it is essential to elucidate the binding mechanism of DNA binding domain (DBD) of transcription factors (TFs), and predict the location of transcription factor binding sites (TFBSs). For an exhaustive search of TFBSs, we have investigated four typical TFs with diverse origins, such as WRKY, PU.1, GLUCOCORTICOID RECEPTOR (GR), and MYC2 by using a newly developed method, KaScape. During KaScape experiments, we identified short sequences (3-4 bases) or "anchoring element" (AE) for the four TFs that dominated the bound population of DNA-DBD binding. We further developed the AEEscape (AE Energy landscape) algorithm to detect and confirm the AE and derived its binding energy landscape for all possible sequences. Our analysis of the energy landscape revealed an energetic funnel around the TFBS, which is related to the AE density gradient in the region surrounding the TFBS. Our results provide novel insights into the mechanism of TF binding to TFBSs.

The current COVID-19 pandemic urges in-depth investigation into proteins encoded with coronavirus (CoV), especially conserved CoV replicases. The nsp13 of highly pathogenic MERS-CoV, SARS-CoV-2, and SARS-CoV exhibit the most conserved CoV replicases. Using single-molecule FRET, we observed that MERS-CoV nsp13 unwound DNA in discrete steps of approximately 9 bp when ATP was used. If another NTP was used, then the steps were only 4 to 5 bp. In dwell time analysis, we detected 3 or 4 hidden steps in each unwinding process, which indicated the hydrolysis of 3 or 4 dTTP. Based on crystallographic and biochemical studies of CoV nsp13 helicases, we modeled an unwinding mechanism similar to the spring-loaded mechanism of HCV NS3 helicase, although our model proposes that flexible 1B and stalk domains, by allowing a lag greater than 4 bp during unwinding, cause the accumulated tension on the nsp13-DNA complex. The hinge region between two RecA-like domains in SARS-CoV-2 nsp13 is intrinsically more flexible than in MERS-CoV nsp13 due to the difference of a single amino acid, which causes the former to induce significantly greater NTP hydrolysis. Our findings thus establish a blueprint for determining the unwinding mechanism of a unique helicase family. O_LIWhen dTTP was used as the energy source, 4 hidden steps in each individual unwinding step after 3 - 4 NTP hydrolysis were observed. C_LIO_LIAn unwinding model of MERS-CoV-nsp13 which is similar to the spring-loaded mechanism of HCV NS3 helicase, except the accumulation of tension on nsp13/DNA complex is caused by the flexible 1B and stalk domains that allow a lag of 4-bp in unwinding. C_LIO_LIComparing to MERS-CoV nsp13, the hinge region between two RecA-like domains in SARS-CoV-2 nsp13 is intrinsically more flexible due to a single amino acid difference, which contributes to the significantly higher NTP hydrolysis by SARS-CoV-2 nsp13. C_LI

Microbial rhodopsin, an important photoreceptor protein, has been widely used in several fields, such as optogenetics, biotechnology, and biodevices etc. However, current microbial rhodopsins are all transmembrane proteins, which both complicates the investigation on the photoreaction mechanism and limits their further applications. Therefore, a suitable mimic for microbial rhodopsin can not only provide a better model for understanding the mechanism, but also can extend the applications. The human protein CRABPII turns out to be a good template for design mimics on rhodopsin, due to the convenience in synthesis and the stability after mutations. Recently, Geiger et al. designed a new CRABPII-based mimic M1-L121E on microbial rhodopsin with the correct 13-cis (13C) isomerization after irritation. However, it still remains a question how similar it is compared with the natural microbial rhodopsin, in particular in the aspect of the photoreaction dynamics. In this article, we investigated the excited-state dynamics of this mimic by measuring its transient absorption spectra. Our results reveal that there are two components in the solution of mimic M1-L121E at PH=8, known as protonated Schiff base (PSB) and unprotonated Schiff base (USB) states. In both states, the photoreaction process from 13-cis (13C) to all-trans (AT) is faster than that from the inverse direction. In addition, the photoreaction process in PSB state is faster than that in the USB state. In the end, we compared the isomerization time of the PSB state with the properties of the microbial rhodopsin, and confirmed that the mimic M1-L121E indeed captures the main feature of the rhodopsin and is a good model of microbial rhodopsin in the photoreaction dynamics. However, our results also reveal significant differences in the excited-state dynamics of the mimic relative to the natural microbial rhodopsin, including the slower PSB isomerization rates in both 13C-AT and AT-13C directions, as well as the unusual USB photoreaction dynamics at PH=8. Such unique properties have not been observed in the natural rhodopsin, which could further deepen the understanding in photoreaction mechanism of the photosensitive proteins.

Avena Sativa phototropin 1 Light-oxygen-voltage 2 domain (AsLOV2) is the model protein of Per-Arnt-Sim (PAS) superfamily, characterized by conformational changes in response to external environmental stimuli. This conformational change is initiated by the unfolding of the N-terminal helix in the dark state followed by the unfolding of the C-terminal helix. The light state is defined by the unfolded termini and the subsequent modifications in hydrogen bond patterns. In this photoreceptor, {beta}-sheets have been identified as crucial components for mediating allosteric signal transmission between the two termini. In this study, we combined microsecond all-atm molecular dynamics simulations and Markov state modeling of conformational states to quantify molecular basis of mutation-induced allostery in the AsLOV2 protein. Through a combination of computational investigations, we determine that the H{beta} and I{beta} strands are the most critical structural elements involved in the allosteric mechanism. To elucidate the role of these {beta}-sheets, we introduced 13 distinct mutations (F490L, N492A, L493A, F494L, H495L, L496F, Q497A, R500A, F509L, Q513A, L514A, D515V, and T517V) and conducted comprehensive simulation analysis. The results highlighted the role of two hydrogen bond Asn482-Leu453 and Gln479-Val520 in the observed distinct behaviors of L493A, L496F, Q497A, and D515V mutants. The comprehensive atomistic-level analysis of the conformational landscapes revealed the critical functional role of {beta}-sheet segments in the transmission of the allosteric signal upon the photoactivation of the light state.

In the traditional scoring functions for protein-ligand binding affinity prediction, the energies of the electrostatic and van der Waals interactions were evaluated (or restricted) by the mathematical expressions of [Formula] and [Formula], respectively. In comparison, the power exponents of distance-based variables as adopted in the present study are not restricted as those in traditional energy terms for atomic interactions. The distance-based variables were integrated using linear regression and artificial neural network to predict the protein-ligand binding affinity or binding energy. The training of the linear, neural network and mixed models was based on the newest data in PDBbind, i.e., PDBbind (v.2024). Estimated according to Pearsons correlation coefficient (R), the best performances of the linear models are 0.700 < R [&le;] 0.800 with the high-quality affinity data, and those of the neural network-based mixed models are 0.800 [&le;] R < 0.900 with the same data. The predictive powers of the best models developed in this study are superior to the sophisticated linear and machine learning-based scoring functions developed before. The results suggest that the distance-based variables with appropriate power exponents may have the ability to improve the prediction of protein-ligand binding affinity with high accuracy. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=143 HEIGHT=200 SRC="FIGDIR/small/680851v1_ufig1.gif" ALT="Figure 1"> View larger version (50K): org.highwire.dtl.DTLVardef@10e86e4org.highwire.dtl.DTLVardef@b9ec82org.highwire.dtl.DTLVardef@565a90org.highwire.dtl.DTLVardef@153c4ef_HPS_FORMAT_FIGEXP M_FIG C_FIG HIGHLIGHTSO_LIBy using the newest data in PDBbind (v.2024) to train the linear, neural network and mixed models, the quantitative distance-energy relationships are further explored and improved to predict the binding affinity of protein-ligand complexes. C_LIO_LIThe power exponents of distance in the traditional energy terms are expanded to characterize the distance-energy relationships accurately at atom level for protein-ligand interactions. C_LIO_LIThe best models are superior to the sophisticated machine learning-based scoring functions developed before. C_LI

Aromatic side-chains (phenylalanine and tyrosine) of a protein flip by 180{degrees} around the C{beta} - C{gamma} axis ({chi}2 dihedral of side-chain) producing two symmetry-equivalent states. The ring-flip dynamics act as an NMR probe to understand local conformational fluctuations. Ring-flips are categorized as slow (ms onwards) or fast (ns to near ms) based on timescales accessible to NMR experiments. In this study, we investigated the ability of the infrequent metadynamics approach to discriminate between slow and fast ring-flips for eight individual aromatic side-chains (F4, Y10, Y21, F22, Y23, F33, Y35, F45) of basic pancreatic trypsin inhibitor (BPTI). Well-tempered metadynamics simulations were performed to observe ring-flipping free energy surfaces for all eight aromatic residues. The results indicate that{chi} 2 as a standalone collective variable (CV) is not sufficient to classify fast and slow ring-flips. Most of the residues needed{chi} 1 (N - C) as a complementary CV, indicating the importance of librational motions in ring-flips. Multiple pathways and mechanisms were observed for residues F4, Y10, and F22. Recrossing events are observed for residues F22 and F33, indicating a possible role of friction effects in the ring-flipping. The results demonstrate the successful application of the metadynamics based approach to estimate ring-flip rates of aromatic residues in BPTI and identify certain limitations of the approach.

Adenosine triphosphate (ATP) is often referred as the energy currency of the cell. Yet, non-invasive, real-time, and quantitative measurement of its concentration in living mammalian cells has been difficult. Here we report an improved fluorescent ATP indicator protein, QUEEN-37C, which is optimized for measuring ATP concentration in living mammalian cells. Absolute value of the ATP concentration can be estimated from the ratiometric fluorescence imaging, and its accuracy was verified by the luciferase assay. Since QUEEN-37C enables the single-cell measurement of ATP concentration, we can not only measure its mean but its distribution in the cell population, which revealed that the ATP concentration is tightly regulated in most cells. We also noted the positive correlations in the ATP concentration among adjacent cells in epithelial cell sheet and mouse embryonic stem cell colonies. Thus, QUEEN-37C would serve as a new tool for the investigation of the single cell heterogeneity of metabolic states.

The clock protein of cyanobacteria KaiC forms a homohexamer with two ring-shaped domains, C1 and C2. These domains undergo several domain-specific conformational transitions and allosterically communicate to generate a circadian rhythm. Interestingly, experiments show a possibility that C2 is independent of C1. However, detailed interplay among them remains elusive. Here we propose a mathematical model, which explicitly considers the interplay. The allostery in KaiC is here modeled to be unidirectional from C2 to C1. We demonstrate that the unidirectional allostery is sufficient for the circadian rhythm by showing the quantitative reproducibility of various experimental data, including temperature dependence of both phosphorylation oscillation and ATPase activity. Based on the present model, we further discuss possible functional roles of the unidirectional allostery particularly in the period robustness against both protein concentration and temperature.

Some calcium channels selectively permeate Ca2+, despite the high concentration of monovalent ions in the surrounding environment, which is essential for many physiological processes. Without atomistic and dynamical ion permeation details, the underlying mechanism of Ca2+ selectivity has long been an intensively studied, yet controversial, topic. This study takes advantage of the homologous Ca2+-selective TRPV6 and non-selective TRPV1 and utilizes the recently solved open-state structures and a newly developed multi-site calcium model to investigate the ion binding and permeation features in TRPV channels by molecular dynamics simulations. Our results revealed that the open-state TRPV6 and TRPV1 show distinct ion-binding patterns in the selectivity filter, which lead to different ion permeation features. Two Ca2+ ions simultaneously bind to the selectivity filter of TRPV6 compared with only one Ca2+ in case of TRPV1. Multiple Ca2+ binding at the selectivity filter of TRPV6 permeated in a concerted manner, which could efficiently block the permeation of Na+. Cations of various valences differentiate between the binding sites at the entrance of the selectivity filter in TRPV6. Ca2+ preferentially binds to the central site with a higher probability of permeation, repelling Na+ to a peripheral site. Therefore, we believe that ion binding competition at the selectivity filter of calcium channels, including the binding strength and number of binding sites, determines Ca2+ selectivity under physiological conditions. Additionally, our results showed that pore helix flexibility and the cytosolic domain of TRPV channels regulate ion permeability.

To elucidate computationally a binding mechanism of a middle-sized flexible molecule, bosentan, to a GPCR protein, human endothelin receptor type B (hETB), a GA-guided multidimensional virtual-system coupled molecular dynamics (GA-mD-VcMD) simulation was performed. This method is one of generalized ensemble methods and produces a free-energy landscape of the ligand-receptor binding by searching large-scale motions accompanied with stably keeping the fragile cell-membrane structure. All molecular components (bosentan, hETB, membrane, and solvent) were represented with an all-atom model, and sampling was carried out from conformations where bosentan was distant from the binding site in the hETBs binding pocket. The deepest basin in the resultant free-energy landscape was assigned to the native-like complex conformation. The obtained binding mechanism is as follows. First, bosentan fluctuating randomly in solution is captured by a tip region of the flexible N-terminal tail of hETB via nonspecific attractive interactions (fly-casting). Bosentan then occasionally slides from the tip to root of the N-terminal tail (ligand-sliding). In this sliding, bosentan passes the gate of the binding pocket from outside to inside of the pocket with accompanying a quick reduction of the molecular orientational variety of bosentan (orientational selection). Last, in the pocket, ligand-receptor attractive native contacts are formed, and eventually the native-like complex is completed. The bosentan-captured conformations by the tip- and root-regions of the N-terminal tail correspond to two basins in the free-energy landscape, and the ligand-sliding corresponds to overcoming a free-energy barrier between the basins.

Circadian rhythm by Cyanobacteria is one of the simplest biological clocks: the clock consists of only three proteins, KaiA, KaiB and KaiC. Their oligomers, KaiA dimer (A2), KaiB tetramer (B4) and KaiC hexamer (C6) oscillate an association- disassociation cycle with 24hr period. In a widely accepted model, the oscillation process is as follows. From the viewpoint of a base unit (C6), C6 homo-oligomer [-&gt;] A2C6 complex [-&gt;] B6C6 complex [-&gt;] AnB6C6 complex (n[&le;]12) [-&gt;]C6 homo-oligomer. In this study, Small-Angle X-ray Scattering, Contrast Matching-Small-Angle Neutron Scattering, Analytical Ultracentrifuge and phosphorylation-analysis PAGE measurements were performed to reveal the kinetics not only of KaiC hexamer but also of all components in a working Kai clock. The complementary analysis disclosed that the oscillation is not the single process as the widely accepted model but composed with synchronized multiple association-dissociation reactions between components. Namely, there are various reactions between components, which proceed simultaneously, in a working Kai-clock.

Glioma is lethal malignant brain cancers, many reports have shown that abnormalities in the behavior of water and ion channels play an important role in regulating glioma proliferation, migration, apoptosis and differentiation. Recently, new studies have suggested that some long noncoding RNAs (lncRNAs) containing small open reading frames (smORFs), can encode small peptides and form oligomers for water or ion regulation. However, because these peptides are difficult to identify, their functional mechanisms are far from being clearly understood. In this study, we used bioinformatic methods and softwares to identify and evaluate lncRNAs in gliomas that may encode small transmembrane peptides. Combining ab initio homology modeling, molecular dynamics simulations and energetic calculations, we constructed a predictive model and predicted the oligomer channel activity of peptides by identifying the lncRNA ORFs. We found that one key hub lncRNA, DLEU1, which contains two smORFs (ORF1 and ORF8) could encode small peptides that form pentameric channels. The mechanics of water and ion (Na + and Cl-) transport through this pentameric channel were simulated. The potential of mean force (PMF) of the H2O molecules along the two ORF-encoded peptide channels indicated that the energy barrier was different between ORF1 and ORF8. The ORF1-encoded peptide pentamer acted as a self-assembled water channel but not as an ion channel, and the ORF8 neither permeating ions nor water. This work provides new methods and theoretical support for further elucidation of the function of lncRNA-encoded small peptides and their role in cancer. Additionally, it provides a theoretical basis for drug development.

General anesthetics are indispensable in modern medicine because they induce a reversible loss of consciousness and sensation in humans. On the other hand, their molecular mechanisms of action have not yet been elucidated. Several studies have identified the main targets of some general anesthetics. The structures of {gamma}-aminobutyric acid A (GABAA) receptors with the intravenous anesthetics such as propofol and etomidate have recently been determined. Although these anesthetic-binding structures provide essential insights into the mechanism of action of anesthetics, the detailed molecular mechanism of how the anesthetic binding affects the Cl- permeability of GABAA receptors remains to be elucidated. In this study, we performed coarse-grained molecular dynamics simulations for GABAA receptors and analyzed the resulting simulation trajectories to investigate the effects of anesthetic binding on the motion of GABAA receptors. The results showed large structural fluctuations in GABAA receptors, correlations of motion between the amino-acid residues, large amplitude motion, and autocorrelated slow motion, which were obtained by advanced statistical analyses. In addition, comparison of the resulting trajectories in the presence or absence of the anesthetic molecules revealed a characteristic pore motion related to the gate-opening motion of GABAA receptors.

The highly infectious SARS-CoV-2 variant B.1.617 with double mutations E484Q and L452R in the receptor binding domain (RBD) of SARS-CoV-2s spike protein is worrisome. Demonstrated in crystal structures, the residues 452 and 484 in RBD are not in direct contact with interfacial residues in the angiotensin converting enzyme 2 (ACE2). This suggests that albeit there are some possibly nonlocal effects, the E484Q and L452R mutations might not significantly affect RBDs binding with ACE2, which is an important step for viral entry into host cells. Thus, without the known molecular mechanism, these two successful mutations (from the point of view of SARS-CoV-2) can be hypothesized to evade human antibodies. Using in silico all-atom molecular dynamics (MD) simulation as well as deep learning (DL) approaches, here we show that these two mutations significantly reduce the binding affinity between RBD and the antibody LY-CoV555 (also named as Bamlanivimab) that was proven to be efficacious for neutralizing the wide-type SARS-CoV-2. With the revealed molecular mechanism on how L452R and E484K evade LY-CoV555, we expect that more specific therapeutic antibodies can be accordingly designed and/or a precision mixing of antibodies can be achieved in a cocktail treatment for patients infected with the variant B.1.617. O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Biomolecular condensation is involved in various cellular processes both functional and dysfunctional. Regulation of the condensation is thus crucial to avoid pathological protein aggregation and to maintain stable cellular environments. Recently, a class of highly charged intrinsically disordered proteins (IDPs), which are called the heat-resistant obscure (Hero) proteins, are shown to protect other client proteins from pathological aggregation. Besides the potential importance of this function, molecular mechanisms for how Hero proteins can protect other proteins from aggregation are not still known. Here we perform multiscale molecular dynamics (MD) simulations of Hero11, one of the Hero proteins, and the C-terminal region of TDP-43, as a target protein of Hero11, at various conditions to examine how they interact with each other. Based on the simulation results, three possible mechanisms have been proposed: (i) TDP-43 and Hero11 in dense phase reduces contacts with each other and shows faster diffusion due to the repulsive Hero11-Hero11 interactions, (ii) the amount of TDP-43 in dilute phase increases and their sizes become greater upon the attractive Hero11-TDP-43 interactions, and (iii) Hero-11 on the surface of small TDP-43 condensates avoids their fusions with the repulsive interactions. We also examine possible Hero-11 structures in atomistic and coarse-grained MD simulations and found disordered Hero-11 tend to assemble on the surface of the condensates, avoiding the droplet fusion effectively. The proposed mechanisms give us new insight into the regulation of biomolecular condensation in the cells and other conditions.

The coronavirus spike protein, which binds to the same human receptor in both SARS-CoV-1 and 2, has been implied to be a potential source of their differential transmissibility. However, the mechanistic details of spike protein binding to its human receptor remain elusive at the molecular level. Here, we have used an extensive set of unbiased and biased microsecond-level all-atom molecular dynamics (MD) simulations of SARS-CoV-1 and 2 spike proteins to determine the differential dynamic behavior of prefusion spike protein structure in the two viruses. Our results indicate that the active form of the SARS-CoV-2 spike protein is more stable than that of SARS-CoV-1 and the energy barrier associated with the activation is higher in SARS-CoV-2. Our results also suggest that not only the receptor binding domain (RBD) but also other domains such as the N-terminal domain (NTD) could play a role in the differential binding behavior of SARS-CoV-1 and 2 spike proteins.

The biological activity of tRNA is closely related to its mechanical folding properties. Although previous studies focused on the folding and unfolding mechanism of tRNA, its kinetics are largely unknown. In this study, combining optical tweezers and molecule dynamics simulations, we characterized the mechanical folding and unfolding processes of a single unmodified Saccharomyces cerevisiae tRNAphe. We identified the intermediates and pathways for tRNA mechanical folding and unfolding in the presence of Mg2+, discovering that the folding/unfolding kinetics of D stem-loop and T stem-loop but not the anti-codon stem-loop significantly affected by their upstream and downstream structures. The cooperative unfolding of the tRNA in the presence of Mg2+ lead to a large hysteresis between the folding and unfolding pathway, and such hysteresis and unfolding cooperativity are significantly reduced by lowering the Mg2+ concentration or mutating the nucleotides forming the elbow structure. Moreover, both steered molecular dynamics simulation and optical tweezers experiment results support that, formation of tertiary interactions in the elbow region increases energy barriers of the mechanical unfolding pathway, including those in between intermediates, and determines the overall unfolding cooperativity. Our studies may shed light on the detailed tRNA chaperone mechanism of TruB and TrmA.

Plexin-semaphorin signaling regulates key processes such as cell migration, neuronal development, angiogenesis, and immune responses. Plexins stand out because they can directly bind with both Rho- and Ras-family small GTPases through their intracellular domains when these GTPases are in their active, GTP-bound states. This binding occurs via intracellular regions which include a Rho-GTPase Binding Domain (RBD) and a GTPase Activating Protein (GAP) segment. Studies have shown that Rho and Ras GTPases play vital roles in plexin signaling and activation. However, the structural dynamics of plexins and GTPases and how these conformational changes affect interactions when plexin is bound with both Ras and Rho-GTPases or bound to only one specific GTPase has remained unclear. In this study, we conducted molecular dynamics (MD) simulations on six distinct plexin-GTPase bound systems to investigate the differences in conformations and dynamics between Plexin-B1 and three GTPases: Rap1b, Rnd1, and Rac1. Our analysis revealed that dynamics with Rac1 are more altered, compared to Rnd1 depending on whether plexins GAP domain is bound or unbound to Rap1b. In addition, we further investigated alterations in network centralities and compared the network dynamics of the Plexin-GTPases complexes, focusing on the differences when Plexin is bound to both Ras (Rap1b) and Rho-GTPases (Rnd1/Rac1) versus when it is bound to only one GTPase. Our study revealed that Rnd1 exhibits stronger and more stable interactions with Plexin-B1 in the absence of Rap1b, while Rac1 shows fewer and less stable connections in comparison. These computational models have features that broadly agree with experimental results from hydrogen-deuterium exchange detected by mass spectrometry (HDX-MS). Such insights provide a better understanding of the molecular mechanisms underlying Plexin-GTPase interactions and the complexities of signaling mechanisms involving GTPases in general.

Diseases spread by mosquitoes lead to death of 700,000 people each year. The main way to reduce transmission is vector control by biting prevention with chemicals. However, the most commonly used insecticides lose efficacy due to the growing resistance. Voltage-gated sodium channels (VGSCs), membrane proteins responsible for the depolarizing phase of an action potential, are targeted by a broad range of neurotoxins, including pyrethroids and sodium channel blocker insecticides (SCBIs). Reduced sensitivity of the target protein due to the point mutations threatened malaria control with pyrethroids. Although SCBIs - indoxacarb (a pre-insecticide bioactivated to DCJW in insects) and metaflumizone - are used in agriculture only, they emerge as promising candidates in mosquito control. Therefore, a thorough understanding of molecular mechanisms of SCBIs action is urgently needed to break the resistance and stop disease transmission. In this study, by performing an extensive combination of equilibrium and enhanced sampling molecular dynamics simulations (3.2 s in total), we found the DIII-DIV fenestration to be the most probable entry route of DCJW to the central cavity of mosquito VGSC. Our study revealed that F1852 is crucial in limiting SCBI access to their binding site. Result explain the role of the F1852T mutation found in resistant insects and the increased toxicity of DCJW compared to its bulkier parent compound, indoxacarb. We also delineated residues that contribute to both SCBIs and non-ester pyrethroid etofenprox binding and thus could be involved in the target site cross-resistance. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=79 SRC="FIGDIR/small/534712v1_ufig1.gif" ALT="Figure 1"> View larger version (27K): org.highwire.dtl.DTLVardef@7912b1org.highwire.dtl.DTLVardef@ba624eorg.highwire.dtl.DTLVardef@8c28b3org.highwire.dtl.DTLVardef@1c42f3c_HPS_FORMAT_FIGEXP M_FIG C_FIG

DNA replication is a high-fidelity information-copying processes which is realized by DNA polymerase (DNAP). The high fidelity was explained on the basis of the well-known kinetic-proofreading mechanism (KPR), under which the so-called fidelity-speed trade-off was studied theoretically. However, numerous biochemical experiments have shown that the high fidelity of DNA replication is achieved due to the initial discrimination of polymerase domain of DNAP, as well as the proofreading of the exonuclease domain of DNAP. This exonuclease-proofreading mechanism (EPR) is totally different from KPR. So the trade-off issues are worth being re-examined under EPR. In this paper, we use the first-passage method recently proposed by us to discuss the possible trade-offs in DNA replication under EPR. We show that there could be no fidelity-speed trade-off under EPR, i.e., the fidelity and the speed can be simultaneously enhanced by EPR in a large range of kinetic parameters. This provides a new perspective to understand the experimental data of the exonuclease activity of T7 DNAP and T4 DNAP. We also show that there exists the fidelity-proofreading cost trade-off, i.e., the fidelity is enhanced at the cost of increasing the futile hydrolysis of dNTP. A possible way to avoid this trade-off is to regulate the rate of DNAP translocation: slowing down the forward translocation (in the presence of the terminal mismatch) can enhance the fidelity without changing the speed and the proofreading cost. Our theoretical analysis offers deeper insights on the kinetics-function relation of DNAP. PACS numbers: 82.39.-k, 87.15.Rn, 87.16.A-